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Numerical  experiments  in  two-  and  three-dimensions,  performed  with  a  high-fidelity  algorithm,  are  em¬ 
ployed  to  explore  plasma-actuator-based  control  of  flow  past  a  stalled  NACA  0015  wing  section.  The  time- 
accurate  force  field  from  a  self-consistent  asymmetric-electrode  arrangement  at  5 kHz  (continuous  wave)  is 
coupled  to  the  Navier-Stokes  equations  at  a  nominal  Reynolds  number  of  45, 000  and  angle-of-attack  of  15°. 
When  the  actuator  is  placed  slightly  downstream  of  the  separation  point,  the  mean  and  asymptotic  flow  re¬ 
sponse  include  formation  of  a  wall-jet  and  inhibition  of  separation.  Features  of  the  unsteadiness  directly  as¬ 
sociated  with  the  forcing,  as  well  as  that  due  to  shedding  occurring  near  the  trailing  edge  are  delineated.  The 
results  suggest  that  both  streamwise  and  normal  (towards  the  wall)  components  of  force  play  important  roles. 
A  simpler  phenomenological  model  is  also  employed  to  examine  flow  control  in  the  full  3-D  setting.  At  higher 
Reynolds  number,  90,  000,  the  flow  becomes  turbulent  immediately  downstream  of  the  leading  edge  and  re¬ 
mains  attached  in  the  mean.  The  effect  of  the  actuator  in  this  situation  is  to  laminarize  the  flow  together  with 
the  formation  of  a  wall-jet.  When  the  actuator  location  is  moved  to  the  point  of  maximum  thickness,  laminar 
separation  is  induced  near  the  leading  edge  with  subsequent  breakdown  to  turbulence.  When  stall  is  observed, 
at  18  deg  angle  of  attack,  placing  the  actuator  near  the  leading  edge  yields  the  anticipated  laminar  wall  jet 
on  the  upper  surface.  The  effect  of  actuator  force  strength  on  separation  suppression  is  examined  in  terms  of 
the  combined  effect  of  transition  to  turbulence  and  near  wall  momentum  enhancement.  Finally,  the  effect  of 
pulsing  is  documented  at  Reynolds  number  45,  000  with  20%  duty  cycle  and  interpulse  periods  of  0.7  and  0.25 
characteristic  times  respectively. 


I.  Introduction 

Plasma-based  techniques  exploiting  electromagnetic  forces  for  flow  control  are  currently  of  considerable  interest. 
Particularly  attractive  properties  are  the  absence  of  moving  parts  and  lack  of  mechanical  inertia,  the  latter  facilitating 
near  instantaneous  deployment  over  a  broad  range  of  frequencies.  Recent  efforts  have  explored  many  possibilities  both 
in  high-speed  as  well  as  low-speed  applications.  The  focus  of  this  paper  is  on  the  latter  speed  regime,  where  striking 
experimental  observations  have  been  obtained  at  atmospheric  pressures  with  glow  discharges1  or  dielectric  barrier 
discharge  (DBD)  based  devices2  operating  in  the  low  radio  frequency  range  of  1  —  10kHz  and  voltage  amplitudes 
of  5  —  10 kV.  References  3-9  include  some  recent  observations  on  separation  control  over  various  airfoil  and  turbine 
sections  under  static  and  oscillating  conditions.  An  overview  of  some  basic  concepts  behind  this  control  technique, 
including  impact  of  duty  cycle,  may  be  found  in  Ref.  10. 

Despite  this  recent  progress,  a  clear  understanding  of  the  physics  of  operation  is  presently  lacking.  A  framework 
for  analysis  was  proposed  in  Ref.  2  by  associating  the  phenomenology  to  the  characteristics  of  the  DBD,  thus  bringing 
to  bear  kinetic  theories  developed  in  plasma  physics  disciplines  in  non-aerospace  contexts.  Further  diagnostic  efforts 
have  been  reported  in  Ref.  11,  which,  based  on  acoustic  measurements,  concluded  that  compressibility  effects  may 
play  a  role  in  momentum  coupling,  though  as  noted  in  Ref.  10  there  is  no  significant  heating  of  the  air. 

Numerical  modeling  of  high-Reynolds  number  flows  with  purely  first-principles  based  approaches  is  a  daunting 
endeavor  because  quantitative  aspects  of  many  of  the  fundamental  processes  remain  either  unknown  or  cannot  be 
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resolved  for  the  parameters  of  interest  because  of  computational  complexity,  especially  for  air  Additionally,  situations 
in  which  flow  control  is  essential,  such  as  stalled  airfoils,  are  characterized  by  flow  transition  to  turbulence,  which 
require  massive  simulation  resources  even  without  control.  These  difficulties  have  fostered  the  development  of  a 
wide  spectrum  of  models  introducing  varying  degrees  of  empirical  elements  into  the  procedure.  Among  the  simpler 
methods  focused  specifically  on  discharge-fluid  coupling,  are  those  of  Roth12  and  Shyy  et  al.13  Roth  associated  transfer 
of  momentum  from  charged  to  neutral  particles  as  effectively  based  on  a  gradient  of  electric  pressure,  which  varies  as 
the  square  of  the  electric  field  magnitude.  A  more  sophisticated  model,  suitable  for  coupling  with  the  fluid  response 
was  proposed  by  Shyy  et  a/,13  who  specify  a  spatially  variable  force  field  suited  for  coupling  to  high-fidelity  numerical 
procedures.  Other  efforts  include  those  of  Corke  et  al 4  where  the  wall-jet  characteristics  are  described  and  compared 
with  experiment  and  of  Voikhov  et  al ,14  who  describe  several  2-D  simulations  with  a  consistent  approach  satisfying 
Gauss’  law. 

The  quest  for  more  sophisticated  tools  utilizing  plasma  kinetic  processes  for  first-principles  simulations  is  being 
pursued  on  several  fronts.  The  problem  requires  self-consistent  solution  of  multi-dimensional  multi-fluid  equations, 
which  implicitly  couple  the  Maxwell  and  Navier-Stokes  equations  and  incorporate  various  inter-molecular  and  elec¬ 
tronic  phenomena.  Recent  efforts  with  varying  levels  of  sophistication  are  reported  in  Refs.  15-18.  Of  interest  in  the 
present  paper  is  the  procedure  developed  by  Roy  et  al  initially  for  1-D  phenomena19  with  subsequent  extensions  to 
2-D  configurations.20-21  In  a  recent  effort.  Ref.  22,  parametric  studies  are  described  examining  the  effect  of  amplitude 
and  shape  of  excitation,  dielectric  constants,  initial  ionization  level,  electrode  shape  and  sequential  ganging  of  actua¬ 
tors.  Key  attributes  of  the  method  are  simultaneous  treatment  of  events  in  both  the  fluid  as  well  as  the  dielectric,  an 
integrated  approach  to  coupling  the  Poisson  equation  and  innovative  methods  of  specifying  boundary  conditions. 

In  the  present  paper,  the  numerical  effort  initiated  in  Ref.  23  to  examine  flow  control  past  a  stalled  NACA  0015 
wing  section  at  a  nominal  Reynolds  number  of  45,  000  and  angle-of-attack  aa  =  15°  (see  Fig.  la)  is  extended  in 
several  different  respects.  First,  the  response  of  the  flow  to  spatially  and  temporally  varying  force  field  obtained  from  a 
first-principles  simulation21-22  is  examined  by  direct  coupling  to  the  2-D  Navier-Stokes  equations.  The  mean  response 
is  qualitatively  compared  to  that  obtained  from  a  previously  developed  phenomenological  model  of  the  average  force 
field.  The  relative  impact  of  the  negative  voltage  phase  of  the  RF  cycle  is  evaluated  by  damping  the  force  during  this 
time  interval.  In  order  to  examine  the  effect  of  Reynolds  number,  angle  of  attack  and  force  magnitude  scaling,  a  3-D 
approach  is  followed  with  a  simplified  force  model.  Finally,  the  effect  of  pulsing  the  force  is  also  explored  by  fixing 
the  duty  cycle  at  20%  at  two  different  interpulse  periods  of  0.7  and  0.25  Tc  respectively,  where  Tc  is  the  characteristic 
time. 


II.  Governing  equations 

The  Navier-Stokes  equations  are  solved  with  body  force  terms.  In  strong-conservation  flux  vector  form  and  general 
curvilinear  coordinates  x  =  x  (£,  y,  Q,y  =  y  (£,  y,  Q,  z  =  z  (£,  y,  £): 

dX  OF,  3G/  dHi  dFv  dGv  dHv  ~ 

-af  +  W  +  lF  +  lK=-sr  +  lF  +  -sr  +  s  (1) 

where,  X  =  {p,  pu,  pv,  pw,  pe}/J  is  the  solution  vector,  J  is  the  transformation  Jacobian  and  F ,  G  and  H  are 
the  flux  vectors  (subscripts  I  and  V  refer  to  inviscid  and  viscous  components  respectively).  All  quantities  are  non- 
dimensionalized  with  reference  values,  length  c  (chord),  velocity  Uref  =  Uoa,  density  pref  =  p^  and  temperature 
Too.  In  terms  of  the  Cartesian  fluxes, 

Fi  =  -j(€xFi  +  £yGi  +  fT/)  Fv  =  —{ixFv  +  £,yGy  +  tizHy)  (2) 

with  similar  expressions  for  the  remaining  flux  vectors.  For  example,  Fj,  Fy  and  S  are: 


pu 

0 

0 

pu 2  +  p 
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where  the  Reynolds  number  is  Re,  Pr  is  the  Prandtl  number  and  M  is  the  Mach  number.  S  =  S/J,  is  the  source 
vector  containing  the  body  force  terms  due  to  the  plasma  actuator.  ec  is  the  electronic  charge  and  the  parameter  Dc 
represents  the  scaling  of  the  electrical  to  inertial  forces: 


Dc 


Pc,ref  &cDref  L ref 
PrefU^ej 


(4) 


Details  of  the  charge  number  density,  pc  and  electric  field  vector  E  =  { Ex ,  Ey ,  Ez }  are  derived  from  models  in¬ 
corporating  various  degrees  of  phenomenological  and  first-principles  components.  In  this  work,  both  approaches  are 
employed  as  described  below,  the  former  in  2-D  to  examine  the  effect  of  the  high-frequency  forcing  on  the  flow  and 
the  latter  in  direct  numerical  simulations  (3-D)  to  examine  effect  with  variation  of  flow  parameters  (Reynolds  number 
and  angle-of-attack)  and  force  pulsing. 


III.  Numerical  Details 


All  simulations  are  performed  at  a  very  low  Mach  number  of  0.1.  Both  2-D  and  3-D  flows  are  simulated,  the 
former  for  computationally  intensive  coupling  with  the  first-principles  unsteady  force  and  the  latter  to  examine  the 
breakdown  process  with  a  phenomenologically  modeled  force.  For  the  3-D  case,  the  spanwise  extent  of  the  domain  is 
taken  to  be  0.2c,  where  c  is  the  chord  length.  A  O-type  mesh  is  employed,  comprised  of  308  x  75  x  145  points  in  the 
x  (streamwise),  y  (spanwise)  and  z  (body-normal)  directions  respectively.  The  grid  is  generated  by  stacking  planes  in 
the  manner  shown  in  Fig.  lb.  A  view  of  each  section  is  shown  in  Fig.  lc,  with  an  enlarged  view  in  Fig.  Id.  The  mesh  is 
stretched  rapidly  in  the  far-field  towards  the  outer  boundary,  which  is  located  30c  from  the  surface  of  the  airfoil.  The 
2-D  grid  is  generated  by  simply  extracting  a  plane  from  the  3-D  mesh.  The  boundary  conditions  are  straightforward. 


a)  Problem  simulated  b)  Overview  of  grid  structure  e)  Sectional  view  d)  Enlarged  view 

Figure  1.  Grid  structure 


No-slip,  zero  body-normal  pressure-gradient  and  isothermal  wall  conditions  are  enforced  at  the  solid-wall.  The  far 
field  is  assumed  to  be  far  enough  away  for  free-stream  conditions  to  be  valid  while  periodic  conditions  are  applied  in 
the  spanwise  direction  as  well  as  at  the  branch  cut  arising  due  to  the  O-type  mesh. 

A  high-order  compact-difference  method  is  employed  to  solve  the  governing  equations.  Derivatives,  <j>',  of  each 
required  quantity,  </>,  are  obtained  in  the  uniformly  discretized  transformed  plane  (0  ij,  Q  with  the  formula: 


Ts(j)'i_l  +</>'  +  Ys(j)'i+1  —  bs 


Pi+ 2 


Pi-2 


Pi+1 


Pi-1 


(5) 


where  Ts,  as  and  bs  determine  the  spatial  properties  of  the  scheme.  All  simulations  described  in  this  work  employ  the 
sixth-order  C6  method,  for  which  Fs  =  1/3,  as  =  14/9  and  bs  =  1/9.  Details  of  the  spectral  characteristics  of  these 
schemes,  and  others  obtainable  from  Eqn.  5,  may  be  found  in  Refs.  24-26. 

Additional  components  are  required  in  the  method  to  enforce  numerical  stability,  which  can  be  compromised  by 
mesh  stretching,  boundary  condition  implementation  and  non-linear  phenomena.  Spurious  frequencies  in  the  solution 
are  systematically  removed  with  a  filter  designed  using  spectral  analysis.  For  any  component  of  the  conserved  vector. 


Otf4>i—l  +  —  ^n=0  c)  [fyi+n 


(6) 


where  is  the  filtered  value.  Relevant  spectral  analyses  may  be  found  in  Ref.  27.  In  the  present  work,  the  8th  order 
filter  with  a/  =  0.2  is  employed  in  the  body-normal  and  streamwise  directions,  while  the  1 0th  order  filter  with 
af  =  0.4  is  chosen  in  the  spanwise  periodic  direction. 
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IV.  Results 


The  baseline  flow  has  been  described  in  a  previous  effort.23  The  key  features  are  reproduced  in  Fig.  2  showing 
iso-levels  of  instantaneous  vorticity  magnitude  colored  by  the  spanwise  component  of  vorticity  (Fig.  2a)  and  instan¬ 
taneous  and  mean  velocity  contours  in  the  mid-span  plane.  The  breakdown  of  the  shear  layer  very  shortly  after 
separation  is  clearly  evident,  generating  smaller  structures  and  spanwise  instabilities.  Although  the  process  is  highly 
three-dimensional,  and  requires  a  direct  numerical  simulation  for  complete  representation,  in  order  to  describe  the 
flow  succinctly,  instantaneous  quantities  will  be  plotted  at  the  mid-span  plane.  Figures  2b  and  c  show  the  instanta¬ 
neous  (/-velocity  and  the  vorticity  magnitude.  The  shear  layer  emanates  from  the  separation  point,  which  occurs  at 
approximately  2%  chord  and  the  flow  may  be  considered  to  be  fully  stalled.  Proceeding  downstream,  the  layer  loses 
its  coherence  as  the  three-dimensional  break-up  process  progresses.  The  flow  beneath  the  shear  layer  is  transitional, 
and  is  characterized  by  very  low  velocities.  Reference  23  contains  additional  description  of  the  no-control  3-D  flow, 
pertinent  aspects  of  which  will  be  reiterated  in  the  context  of  the  new  results  below  where  warranted. 


a)  3-D  structure  -  isolevs 


b)  instantaneous  //-velocity  (mid-span) 

Figure  2.  Structure  of  3-D  flowfield  without  control 


c)  mean  //-velocity 


The  two-dimensional  baseline  flow  is  rather  different  than  in  3-D  because  of  the  absence  of  a  spanwise  breakdown 
mechanism.  Contours  of  select  instantaneous  and  mean  quantities  are  plotted  in  Fig.  3.  In  this  case,  shed  vortices 
maintain  their  coherence  in  an  unphysical  fashion,  giving  rise  to  a  sequence  of  large  vortex  pair  structures  that  are 
shed  downstream.  Peak  instantaneous  velocities.  Fig.  3a,  are  considerably  larger  in  localized  regions,  reaching  almost 
twice  those  observed  in  the  3-D  simulations.  In  a  consistent  fashion,  instantaneous  pressure  contours  exhibit  coherent 
features  in  which  extrema  are  much  higher  than  in  the  full  3-D  simulation.  The  mean  u  contours.  Fig.  3c,  shows 
significant  qualitative  and  quantitative  differences  from  both  instantaneous  2-D  and  mean  3-D  counterparts,  with  a 
lower  incline  of  the  separated  shear  layer  relative  to  the  surface  of  the  airfoil. 


a)  Instantaneous  u  velocity 


b)  Instantaneous  pressure 

Figure  3.  Structure  obtained  with  2-D  simulations  without  control 


c)  Mean  u  velocity 


A.  Coupling  of  force  field  obtained  from  first-principles  approach 

1.  Time-varying  nature  of  force  field 

The  plasma  effect  is  coupled  to  the  Navier-Stokes  equations  through  the  body  force,  which  depends  on  the  net  space 
charge  and  electric  field  distributions.  The  charge  density,  pc  and  the  electric  field  E  are  obtained  in  a  self-consistent 


4  of  20 


American  Institute  of  Aeronautics  and  Astronautics 


fashion  following  the  procedure  described  in  Ref.  21,22.  Briefly,  the  two-dimensional  three-species  collisional  plasma- 
sheath  model  includes  the  charge  and  momentum  continuity  equations,  and  Gauss’  law  for  electric  potential: 


dna 

~dt 


+  — naVaj  =  nez  -  rnen.i 

OXj 


(7) 


Tl0  Vfyj 


=  sgn  {e)napc 


dtp 

dxj 


-Dn 


dna 

dxj 


e 


e  (ne  -  Hi) 


(8) 

(9) 


The  charged  particle  subscript  a  refer  to  electrons,  e  or  ions,  i  respectively,  n  is  the  number  density,  V:l  the  velocity 
vector  and  z  and  r  are  the  ionization  and  recombination  coefficients.  Since  air  chemistry  for  the  pertinent  processes 
in  the  discharge  remain  poorly  understood,  several  assumptions  are  made  to  obtain  representative  features  of  the  time- 
varying  force  field.  First,  the  distribution  and  variation  of  charge  and  electric  field  are  obtained  under  the  assumption 
that  the  working  gas  is  helium  and  secondary  emission  is  currently  ignored.  The  scope  of  this  effort  is  thus  on 
developing  the  coupling  framework  and  examining  fluid  response  to  a  force  fluctuating  at  kHz  frequency.  The  effort 
to  lift  limitations  in  current  thermo-chemical  data  for  the  environment  of  interest  is  being  pursued  separately.  A  degree 
of  quantitative  accuracy  is  reintroduced  through  empiricism  by  scaling  the  force  to  match  peak  amplitudes  obtained 
from  phenomenological  considerations  in  Refs.  13,23.  The  electron  temperature  is  assumed  to  be  uniform  at  leV, 
while  ions  and  neutrals  are  assumed  to  be  in  thermal  equilibrium  at  300/\ .  Electron  diffusion  is  obtained  from  the 
Einstein  relation  while  the  ion  diffusion  coefficient,  mobility,  ionization  and  recombination  rate  coefficients,  which  are 
functions  of  the  electric  field  magnitude  and  the  pressure,  and  taken  from  the  literature.28 

The  above  equations  are  solved  with  the  multiscale  ionized  gas  (MIG)  flow  code,  developed  at  the  Computational 
Plasma  Dynamics  Laboratory  at  Kettering  University.  The  method  is  based  on  a  versatile  finite-element  (FE)  procedure 
adapted  from  fluid  dynamics  to  overcome  the  stiffness  of  the  equations  generated  by  multi-species  charge  separation 
phenomena.  A  2D  bilinear  finite  element  formulation  is  chosen  with  4th  order  Runge-Kutta  time  marching.  The 
solution  process  consists  of  two  steps.  The  first  solves  the  equations  for  charge  and  electric  field  simultaneously.  In  the 
second  step,  the  force  so  obtained  is  transferred  to  the  airfoil  after  rotation  and  scaling  and  constitutes  the  source  term 
in  Eqn.  3,  with  pc  =  ( rii  —  ne).  The  MIG  code  also  solves  for  the  self-consistent  fluid  response,  which  is  not  factored 
in  the  present  analysis  (see  Ref.  22).  This  implicitly  assumes  that  the  near-wall  local  fluid  neutral  velocity  does  not 
influence  the  distribution  of  electric  parameters.  This  requires  that  the  fluid  density  and  pressure,  or  collisionality,  are 
not  much  different  from  those  employed  in  the  plasma  calculation.  This  is  a  reasonable  expectation  for  the  low-speed 
atmospheric  pressure  incompressible  flows  of  interest  here.  Further  details  may  be  found  in  Refs.  22,29-31. 

The  equations  are  solved  in  the  region  surrounding  a  small  flat  surface  which  demarkates  the  ambient  medium 
from  the  dielectric.  The  electrodes  are  placed  in  a  staggered  manner  as  shown  in  Fig.  4a,  with  the  exposed  electrode 
(12mm  long)  being  located  upstream  of  the  embedded  electrode  (12.7 mm  long).  The  vertical  distance  between  the 
electrodes  is  1mm  and  a  0.2 mm  horizontal  overlap  exists  between  the  two.  The  exposed  electrode  is  excited  at  a 
signal  of  5 kHz  and  2 kV  rms  amplitude.  The  dielectric  constant  of  the  fluid  is  assumed  to  be  1.0055eo  while  that  of 
the  solid  is  assumed  to  be  3.5e0,  where  eD  refers  to  the  property  of  free  space.  The  mesh  employed,  shown  in  Fig.  4b, 
consists  of  103  x  106  nodes,  clustered  near  the  electrodes  in  the  vertical  direction,  and  near  electrode  edges  in  the 
streamwise  direction.  The  vertical  height  of  the  simulated  domain  extends  20mm  above  the  surface  and  1  ram  beneath 
it.  Inside  the  dielectric,  the  only  component  of  current  is  associated  with  the  change  of  electric  field  displacement.  At 
the  interface,  this  value  is  matched  with  that  associated  with  electron  flow  in  the  fluid  towards  the  surface  on  which 
charge  accumulation  is  allowed  to  occur.  Details  of  the  variation  of  electromagnetic  fields  at  different  phase  angles  are 
provided  in  Ref.  21.  Charge,  electric  field  and  force  vector  distributions  at  the  peak  of  the  forward  phase  ( tp  =  tt/2) 
are  shown  in  Fig.  4  c  and  d  respectively.  The  charge  distribution  consists  of  regions  of  net  positive  as  well  as  negative 
ranges.  At  the  instant  shown,  the  negative  charge  is  restricted  to  a  very  small  region  at  the  dielectric  surface  above 
the  embedded  electrode.  In  most  regions  away  from  the  surface,  the  ion  concentration  is  larger  and  the  net  charge 
is  positive.  Electric  field  vectors  at  the  same  instant.  Fig.  4d,  are  consistent  with  the  charge  through  Gauss’  law. 
Large  amplitudes  are  evident  near  the  trailing  edge  of  the  exposed  electrode.  At  the  dielectric  surface,  the  direction 
of  the  field  changes  relatively  abruptly  from  upward  (into  the  flow)  to  pointing  downward.  The  magnitude  diminishes 
rapidly  in  the  ^-direction,  together  with  a  reversal  of  direction.  Thus  upstream  pointing  vectors  are  evident  aft  of  the 
exposed  electrode.  The  body  force  field  vectors  are  presented  in  Fig.  4e.  Since  the  charge  is  positive  in  most  regions 
at  this  phase  angle,  the  force  and  electric  field  are  generally  similar.  The  mean  force  field  is  shown  in  Fig.  4f.  In  the 
crucial  region  near  the  trailing  edge  of  the  exposed  electrode,  the  pattern  is  generally  similar  to  that  in  Fig.  4e  with 
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b)  Mesh  employed 
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Figure  4.  Features  of  self-consistent  simulation  of  plasma  actuator  on  flat  plate  configuration 
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a  predominantly  downward/streamwise  orientation.  In  other  regions,  the  force  exhibits  components  away  from  the 
surface  (over  the  exposed  electrode),  as  well  as  towards  it  (above  the  dielectric).  The  impact  of  the  difference  between 
the  variation  of  Fig.  4f  and  the  phenomenological  models  of  Refs.  13, 23,  which  assume  a  triangular  region  of  force, 
is  not  presently  known  for  all  parameters.  It  will  be  shown  later  however,  that  for  the  cases  explored  in  this  study,  the 
average  response  to  the  unsteady  force  bears  great  resemblence  to  that  obtained  with  the  phenomenological  model, 
though  quantitative  differences  exist  for  similarly  scaled  force  fields. 

The  superposition  of  time-varying  force  on  the  airfoil  mesh  is  shown  Fig.  5  at  several  phase  angles  of  the  input 
signal.  The  manner  in  which  this  transformation  was  conducted  is  as  follows.  The  downstream  edge  of  the  exposed 
electrode  was  mapped  to  a  selected  point  on  the  airfoil.  Unless  otherwise  stated,  this  point  lies  slightly  downstream 
of  separation  (XQ  =  x0/c  =  0.028 ,ZQ  =  z0/c  =  0.034)  observed  at  Re  =  45,000,  aa  =  15°.  At  this  point,  a 
local  orthogonal  coordinate  system,  x'  —  z'  is  defined  such  that  x'  is  tangent  to  the  surface.  The  rectangular  domain 
of  Fig.  4b  is  mapped  to  this  local  system  by  simple  rotation  operations.  The  domain  of  active  force  is  then  scaled 
by  a  factor  of  100  so  that  the  length  and  height  of  the  region  are  roughly  3%  and  2%  chord  respectively,  similar  to 
that  employed  with  the  phenomenological  model.23  At  each  point  on  the  airfoil  grid,  the  force  contribution  is  then 
computed  through  bilinear  interpolation,  with  area-weighting  to  ensure  matching  of  integrated  force.  This  process  is 
carried  out  for  each  phase  angle  and  the  force  field  on  the  airfoil  is  stored  for  subsequent  use  as  below.  The  amplitude 
of  the  force  is  first  normalized  by  peak  value  in  the  cycle  -  this  greatly  facilitates  subsequent  magnitude  control  through 
the  coefficient  Dc  defined  above.  Following  Ref.  23,  the  nominal  value  for  Dc  is  2400,  which  is  obtained  by  assuming 
a  charge  density  of  1011/cm3,  peak  electric  field  magnitude  of  226.27 kV/cm  (both  taken  from  Ref.  13),  chord  lenth 
of  12.7cm  and  Reynolds  number  45, 000  at  atmospheric  conditions. 

At  each  time-instant  in  the  Navier-Stokes  simulation,  the  phase  angle  within  a  cycle  is  computed  from  knowledge 
of  the  period  of  the  cycle,  and  the  time  at  which  the  control  was  initiated.  It  is  then  a  simple  matter  to  linearly 
interpolate  between  known  (stored)  adjacent  phase  angles  obtained  from  the  above  spatial  and  temporal  integration 
process.  In  the  cases  described  below,  the  non-dimensional  time-step  was  fixed  at  5  x  1  O  '3.  Since  the  non-dimensional 
period  of  the  5kHz  signal  is  Ta  =  4.7  x  10~3,  about  95  time-steps  span  a  single  cycle. 

The  force  field  variation  over  a  cycle  is  shown  in  Fig.  5.  In  the  “positive”  segment,  i.e.,  where  the  input  voltage 
is  positive,  the  force  field  very  near  the  surface  at  Xa,  Za,  is  directed  downward  towards  the  reference  point.  Slightly 
downstream,  the  force  has  a  significant  streamwise  component  as  well.  Proceeding  outwards  from  the  wall,  the  force 
obtains  an  upward-downstream  orientation.  The  amplitudes  rise  rapidly  with  time  though  a  small  lag  is  observed 
between  magnitude  and  phase-  note  the  relatively  larger  force  field  at  (f>  =  07r/lO  compared  to  that  at  <f>  =  47t/10 
and  the  observation  that  the  force  field  at  f  =  tt  is  non-zero.  On  the  “negative”  stroke,  the  amplitude  of  the  force 
is  relatively  smaller.  The  force  is  predominantly  towards  the  surface,  though  significant  upstream  components  are 
evident  above  the  embedded  electrode. 

2.  Response  to  time-varying  force:  unsteady  flowfield 

When  the  flow  is  subjected  to  the  above  body  force,  the  response  may  be  classified  into  an  initial  response  during 
which  the  large  separated  region  is  eliminated,  and  an  asymptotic  periodic  state.  The  initial  response  is  shown  in  terms 
of  (/-velocity,  pressure  and  spanwise-vorticity  lov,  in  Fig.  6  -  the  flow  time  marked  is  relative  to  the  onset  of  control. 
The  initial  frame,  T  =  0.01TC  shows  the  typical  structure  observed  in  two-dimensional  simulations  of  stall  -  the  flow  is 
characterized  by  large  vortical  structures  that  are  convected  downstream  without  breakdown  because  of  the  absence  of 
a  3-D  instability  mechanisms  in  the  2-D  simulation.  The  control  technique  suppresses  the  formation  of  new  structures 
that  ordinarily  would  have  formed.  The  development  of  the  apparent  attachment  is  relatively  monotonic,  starting  near 
the  point  of  application  of  the  body  force  and  progressing  downstream.  The  initial  signal  reaches  the  trailing  edge  at 
T  ~  1.5TC,  when  the  original  large  vortical  structures  have  convected  downstream  of  the  wing  section.  The  situation 
at  an  instant  after  several  characteristic  times  have  elapsed,  ,  T  =  5.3 Tc,  indicates  an  attached  flow  over  most  of 
the  airfoil,  though  an  unsteady  shear  layer  is  evident  near  the  traililng  edge.  The  acoustic  impact  of  the  actuator  is 
most  obvious  near  the  point  of  application  as  observed  in  the  pressure  contours.  Although  the  acoustic  signal  yields 
large  relative  fluctuations  compared  to  the  no-control  case,  its  amplitude  diminishes  with  distance,  an  effect  which 
is  enhanced  by  the  coarser  mesh  resolution  away  from  the  surface.  Further  details  of  the  unsteadiness  observed  in 
the  pressure  field  are  described  below.  Vorticity  contours  also  demonstrate  the  convection  of  the  initial  structures, 
and,  at  asymptotic  state  indicate  clearly  the  development  of  organized  structures  in  the  shear  layer  separating  near  the 
trailing  edge.  With  this  force  magnitude,  flow  attachment  is  not  complete  in  this  sense.  It  should  be  noted,  following 
the  discussion  in  Ref.  23  and  in  the  context  of  Fig.  3,  that  the  evolution  of  naturally  induced  unsteady  features  such 
as  these  are  not  generally  correctly  reproduced  in  2-D  simulations  because  of  the  absence  of  a  spanwise  breakdown 
mechanism. 
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Figure  5.  Force  field  on  airfoil  at  different  phases  in  a  single  cycle  of  excitation  after  asymptotic  state  is  reached 
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Figure  6.  Initial  transient  and  unsteady  long  term  response  with  spatio-temporally  varying  force 


The  nature  of  the  unsteadiness  encountered  within  an  individual  cycle  are  described  in  Fig.  7.  Surface  static 
pressures  at  the  leading  edge,  mid-chord  and  at  the  trailing  edge  are  plotted  in  Fig  7a  at  intervals  which  are  roughly 
2  actuator  periods  apart.  Large  oscillations  are  evident  near  the  leading  edge,  which,  of  the  three  stations  plotted  is 
closest  to  the  actuator.  The  excursions  from  the  mean  diminish  downstream.  The  fine  scale  nature  of  the  signal  is 
shown  in  Fig.  7b,  where  pressure  is  plotted  at  every  computed  time  instant  for  several  actuator  excitation  periods.  The 
dominant  frequency  near  the  leading  edge  is  evidently  much  higher  than  that  observed  downstream.  Fourier  analysis, 
Fig.  7c,  indicates  that  near  the  leading  edge,  the  primary  frequency  is  that  of  the  actuator,  while  the  smaller  values  at 
the  other  locations  are  associated  with  shedding  of  large  scale  structures,  described  earlier.  The  response  of  velocity 
profiles  in  the  field  at  several  phase  angles  within  one  actuator  cycle  is  exhibited  in  Fig.  7d.  Three  profiles  are  shown, 
just  downstream  of  the  actuator  ( X  =  0.11),  mid-chord  ( X  =  0.5)  and  at  the  trailing  edge  ( X  =  1)  respectively. 
Near  the  actuator,  the  existence  of  the  wall  jet  like  structure,  with  peak  velocities  of  about  2.4Li'oc,  is  prominent. 
Examination  of  the  inset  reveals  that  the  peak  velocity  varies  by  about  4%  within  the  cycle.  At  mid-chord,  the  fullness 
of  the  profile  diminishes,  the  peak  decreases  to  1  .5La,  while  it  location  moves  away  from  the  wall.  At  the  trailing 
edge,  the  velocity  achieves  negative  values  as  the  profile  cuts  across  a  coherent  structure  that  develops  near  the  trailing 
edge  (see  Fig.  6). 


3.  Response  to  time-varying  force:  time-averaged  flowfield 

Aspects  of  the  flow  field  described  above,  averaged  over  several  characteristic  times  after  the  asymptotic  state  is 
reached  are  shown  in  Fig. 8.  Close  examination  of  the  solution  indicates  that  the  point  of  separation  point  moves 
downstream  relative  to  the  no  control  case,  but  only  to  about  x/c  =  0.4,  which  is  beyond  the  maximum  thickness 
location.  However,  the  reversed  flow  layer  remains  restricted  to  the  near  wall  region  for  significant  distance,  and  there 
is  little  displacement  impact  on  the  outer  flow.  Gross  manifestations  of  separation,  including  shear  layer  detachment 
from  the  surface,  become  apparent  only  much  further  downstream,  x/c  ~  0.85,  similar  to  the  experimental  visual¬ 
izations  reported  in  Ref.  12.  The  mean  pressure  contours.  Fig.  8b  are  characteristic  of  an  attached  flow,  with  acoustic 
signals  noted  earlier  in  the  instantaneous  contours  having  been  averaged  out.  The  development  of  the  wall-jet  yields  a 
reversal  of  sign  of  the  sole  vorticity  component  (t oy)  pointing  out  of  the  plane  of  the  figure. 
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4.  Response  to  steady  averaged  force 


Phenomenological  models,  such  as  those  described  in  Refs.  13, 23  employ  a  steady  averaged  force  to  determine  flow 
response.  Even  when  unsteady  force  variation  is  available  as  above,  the  computational  expense,  associated  principally 
with  the  far  smaller  time-step-size  required  to  resolve  each  actuator  period,  makes  it  attractive  to  employ  a  time- 
averaged  force  to  determine  response.  Such  an  approach  does  not  resolve  the  acoustic  phenomena,  but  the  observation 
of  relatively  small  variation  in  velocity  within  the  actuator  cycle,  Fig.7d,  suggests  that  many  key  features  may  be 
accurately  computed  with  an  averaged  force.  For  practical  reasons  therefore,  in  this  section  the  averaged  response  to 
the  unsteady  force  (§IV.  A. 3)  is  compared  with  that  obtained  with  the  steady  average  force,  Fig.5w,  due  to  the  actuator. 

The  asymptotic  response  of  the  flow  to  the  constant  averaged  force  is  unsteady,  because  of  the  manifestation  of 
separation  near  the  trailing  edge  and  subsequent  shedding.  Consequently  the  time-mean  response  is  described  in  Fig. 9 
which  depicts  the  velocity,  pressure  and  vorticity  fields  obtained  by  averaging  the  flow  over  several  characteristic  times 
after  an  asymptotic  state  is  reached.  Comparison  of  each  quantity  in  Fig. 8  and  Fig. 9  indicates  that  the  mean  response 
is  very  similar  suggesting  that  the  detailed  unsteady  characteristics  of  the  body  force  are  not  critical  in  determining 
the  mean  features  of  the  response.  A  more  quantitative  comparison  of  the  two  forcing  approaches  is  shown  in  Fig.  10. 
The  surface  pressure  coefficient.  Fig.  10a,  shows  a  sharp  peak  near  the  point  of  application  of  the  actuator,  but  relaxes 
to  values  typical  of  the  suction  pressure.  Figure  10b  and  c  compare  averaged  u-velocity  profiles  at  several  different 
streamwise  locations.  In  the  immediate  vicinity  of  the  actuator,  X  ~  0.03,  a  small  negative  flow  region  is  evident  very 
near  the  surface,  while  a  wall  jet  with  peak  u-velocity  of  about  3  occurs  in  the  outer  region.  Both  forcing  approaches 
show  the  evolution  of  the  wall  jet,  and  the  subsequent  development  of  a  velocity  deficit  downstream  of  the  trailing  edge 
and  its  subsequent  diffusion  in  the  wake,  where  note  that  the  vertical  scale  has  been  significantly  expanded.  Figure  10 
thus  reiterates  the  similar  mean  response  obtained  from  the  time-accurate  and  averaged  force  fields  inside  the  flow 
and  on  the  surface.  Differences  observed  in  the  velocities  are  modest.  For  example,  at  the  trailing  edge,  the  peak 
discrepancy  is  less  than  5%. 

This  analysis  also  provides  insight  into  the  strengths  and  limitations  of  the  technique  typically  employed  in  phe¬ 
nomenological  approaches.  The  various  models  explored  in  Ref.  23  assume  a  triangular  force  distribution  within 
which  different  orientations  were  chosen.  Fields  with  significant  streamwise  components  were  found  to  be  most  ef¬ 
fective  with  steady  asymptotes  while  body  normal  (towards  the  surface)  components  also  yielded  attached  flows  in  the 
mean.  None  of  the  chosen  distributions  has  the  complex  features  obtained  with  the  first-principles  approach  (Fig. 5). 
Nevertheless,  the  predominant  components  are  streamwise  and  towards  the  surface,  with  relatively  smaller  compo¬ 
nents  upstream  and  away  from  the  surface.  Another  point  of  difference  has  its  basis  in  the  observation  that  for  the 
same  value  of  Dc,  at  most  points  in  the  cycle  the  force  amplitude  is  smaller  than  that  with  the  phenomenological 
approach.  Despite  these  differences,  the  response  is  therefore  relatively  similar,  suggesting  that  while  quantitative 
accuracy  and  acoustic  analysis  demands  first-principles  simulations,  a  broad  preliminary  indication  of  flow  sensitivity 
to  applied  body  force  can  be  conveniently  obtained  through  the  phenomenological  model. 
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Figure  9.  Mean  response  of  flow  to  averaged  force 
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5.  Effect  of  force  during  negative  stroke 

The  relative  impact  of  different  phase  ranges  of  each  individual  cycle  of  the  actuator  excitation  has  been  subject  of 
scrutiny.  Reference  1 1  for  example  divides  the  cycle  into  “negative-going”  and  “positive-going”  segments  depending 
on  the  slope  of  the  voltage.  The  force  field  described  in  Fig. 5  suggests  another  natural  demarkation  into  ranges 
depending  on  the  sign  of  the  voltage.  In  particular,  the  effect  of  the  “negative”  segment  in  the  phase  range  ir  <  <f>  <  2n 
is  of  interest,  since  the  body  force  exhibits  a  significant  upstream  component  in  addition  to  the  body  normal  component. 
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a)  Cp  b)  u  on  upper  surface  C)  u  at  t.e  and  downstream 

Figure  10.  Comparison  of  mean  response  of  flow  to  time-accurate  versus  averaged  force  field 

Results  averaged  after  an  asymptotic  state  is  reached,  are  shown  in  Fig.  1 1  for  force  fields  in  which  the  negative  stroke 
is  damped  to  varying  extents.  In  Fig.  11a,  results  with  the  original  model  are  displayed  while  in  Figs. lib  and  c 
the  reverse  stroke  is  damped  to  50%  and  10%  respectively.  Although  the  effect  of  reducing  the  reverse  stroke  to 
half  its  normal  value  is  not  significant,  the  near  complete  elimination  results  in  a  significant  degradation  of  control 
effectiveness.  In  this  case,  mean  separation  occurs  at  x/c  ~  0.13,  and  the  displacement  effect  of  separation  increases 
sharply.  The  separated  shear  layer  then  exhibits  vortical  structure  formation  similar  to  that  shown  in  Fig. 3.  In  addition 
to  the  upstream  component,  the  negative  stroke  is  also  characterized  by  a  downwards  (toward  the  surface)  component. 
These  results  suggest  that  in  the  absence  of  turbulence,  this  latter  component  during  the  reverse  stroke  is  an  important 
element  of  the  control  effectiveness. 


a)  No  damping 


b)  Negative  stroke  damped  50% 

Figure  11.  Effect  of  damping  negative  stroke  on  mean  ujv 


c)  Negative  stroke  damped  90% 


B.  Effect  of  Reynolds  number  on  control 

All  calculations  described  above,  as  well  as  those  in  Ref.  23,  utilized  a  Reynolds  number  of  45, 000  and  aa  =  15°. 
In  order  to  examine  the  effect  of  Reynolds  number  on  control,  simulations  were  performed  at  Re  =  90,  000  and  aa 
values  of  15°  and  18°  respectively.  In  the  absence  of  control,  the  results  described  below  show  that  the  3-D  turbulence 
structure  is  crucial  in  determining  the  flow  field:  for  example,  at  Re  =  90, 000,  alphaa  =  15°,  the  2-D  simulations 
indicate  laminar  separation  and  formation  of  large  vortical  structures,  while  the  3-D  flow  is  turbulent  and  attached. 
These  calculations  are  thus  performed  in  full  3-D  (direct  numerical  simulations).  However,  when  control  is  applied, 
the  expense  of  the  first-principles  coupling  described  above  becomes  prohibitive  because  of  the  small  time-step-size 
required  to  resolve  each  cycle.  To  alleviate  this  situation,  the  phenomenological  approach  of  Refs.  13, 23  is  employed 
to  obtain  the  body  force  -  this  also  facilitates  comparison  with  Re  =  45,  000  results  of  Ref.  23.  The  results  presented 
earlier  suggest  that  the  key  aspects  of  the  response,  including  the  development  of  the  wall-jet  and  attachment  process, 
are  similar  with  phenomenological  and  first-principles  approaches. 

1.  Flowfield  at  Re  =  90,  000,  aa  =  15°  -  No  control 

No-control  results  at  the  higher  Reynolds  number  are  shown  in  Fig.  12  in  terms  of  instantaneous  velocity  (Fig.  12a), 
vorticity  contours  (Fig.  12b)  and  pressure  (Fig.  12c)  contours  on  the  mid-span  plane.  Mean  velocity  and  vorticity 
contours  are  presented  in  Fig.  13a  and  b,  while  the  mean  surface  pressure  profile  is  compared  with  that  at  Re  =  45,  000 
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in  Fig.  13c.  These  results  indicate  that  the  flow  becomes  turbulent  very  quickly  downstream  of  the  leading  edge,  and 
remains  attached  in  the  mean.  Several  conclusions  are  readily  made  by  comparing  with  Fig. 2,  which  displays  flow 
features  at  Re  =  45, 000.  The  shear  layer  emanating  from  the  separation  point  at  lower  Reynolds  number  is  absent 
as  is  the  large  separated  region  beneath  it.  The  surface  pressure  profile  variation,  Fig.  13c,  shows  the  classic  flat  shape 
at  the  stalled  lower  Reynolds  number  case,  while  the  attached  flow  obtained  at  higher  Reynolds  number  is  yields  a 
characteristic  suction  profile,  with  a  higher  suction  peak. 


b)M 

Figure  12.  Instantaneous  flow  field  features  at  Re  =  90,  000  without  control 


a)  u-velocity  b)  vorticity  c)  Surface  pressure 

Figure  13.  Mean  flow  field  features  at  Re  =  90,  000,  aa  =  15°  without  control 


2.  Brief  description  of  phenomenological  force 

The  effect  of  the  actuator  on  the  flow  was  examined  with  the  phenomenological  force  model  described  in  Ref.  23.  The 
particular  variant  denoted  “Case  2”  is  employed.  In  this,  the  body  force  is  obtained  from 

F  =  Dc6AtapcE5  (10) 

where  9  is  the  frequency  of  the  applied  voltage.  At  is  an  effective  duty  cycle  to  recover  a  mean  force,  a  is  a  factor 
for  collision  efficiency,  pc  is  the  charge  density,  E  is  the  electric  field  vector  and  S  constrains  the  region  of  force 
as  below.  The  charge  density  is  assumed  to  be  uniform  in  the  region  of  interest.  The  variation  in  electric  field  is 
approximated  through  several  assumptions.  The  magnitude  is  assumed  to  vary  linearly  in  x  and  z  measured  from  the 
actuator  position,  \E\  =  Efl  -  k-\ x  —  k%z.  The  peak  value  (Ea  =  V/d,  where  V  is  the  applied  voltage  and  d  is  the 
streamwise  distance  between  the  two  electrodes)  occurs  at  the  origin.  A  further  assumption  is  made  that  the  field  value 
is  at  or  above  breakdown  level  {Ef)  in  the  plasma  region,  so  that  S  =  1  when  \E\  >  E b  and  <5  =  0  otherwise.  By 
choosing  the  streamwise  and  body  normal  extents  ( bd  and  ad  respectively),  the  values  of  ki  =  ( Ea  —  Ef)/bd  and 
=  ( En  —  Ef, ) / a, i  are  determined.  In  the  frame  of  reference  of  the  device,  E  is  then  obtained  from: 

s=F^S"M+K,k'^  <n> 

Here,  £  and  fj  are  unit  vectors  aligned  respectively  along  and  normal  to  the  exposed  electrode,  which  is  assumed  to 
lie  flush  on  the  surface  of  the  body.  For  the  airfoil,  this  coordinate  system  is  located  at  the  tangent  to  the  surface 
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at  the  leading  edge  of  the  electrode.  In  all  cases,  a  single  actuator  is  assumed,  though  its  location  and  strength  are 
varied.  The  body  force  orientation  is  adjusted  by  choosing  k %  and  K.q .  In  this  case,  values  of  1  and  —1  respectively 
were  employed  (vertical  force  component  points  towards  the  surface).  The  values  of  9  =  3000,  At  =  67 ps,  a  =  1, 
pc  =  10 11 /cm3  are  taken  from  Ref.  13.  The  nominal  values  chosen  for  the  parameters  describing  strength  and  extent 
are  Dc  =  2400  (same  as  employed  in  the  results  of  §  IV. A),  ad  =  0.018  are  bd  =  0.024,  where  a  and  b  are  normalized 
by  chord  length.  The  force  field  obtained  when  the  actuator  is  placed  near  the  leading  edge  with  origin  at  Xa  =  0.028, 
Za  =  0.034,  is  shown  in  Fig.  14a  at  the  mid-span  location.  The  force  has  no  gradients  in  the  spanwise  direction. 

3.  Effect  of  actuator  onflow  at  Re  =  90, 000,  aa  =  15° 

Even  though  the  flow  is  attached,  it  is  illustrative  to  examine  the  effect  of  control  because  of  its  implication  from  the 
standpoint  of  laminarization  techniques  as  well  as  off-design  performance.  To  this  end,  the  actuator  was  positioned 
separately  at  two  locations:  the  first  near  the  leading  edge  as  above,  and  the  second  at  mid-chord.  Results  obtained 
with  the  first  actuator  location  are  shown  in  Fig.  14b  through  d.  Over  most  of  the  wing  section,  the  flow  becomes 
attached  and  laminar.  The  wall-jet  is  a  prominent  feature  as  before  -  instantaneous  and  mean  results  are  therefore 
similar  with  the  exception  of  the  region  downstream  of  the  trailing  edge  where  the  shear  layer  becomes  unstable. 
Fig.  14c  compares  (/-velocity  profiles  at  mid-chord  without  and  with  control  for  the  two  Reynolds  numbers  examined. 
Without  control,  the  flow  is  stalled  at  Re  =  45,  000  and  the  mean  w-velocity  is  negative.  At  Re  =  90,  000,  the  velocity 
profile  is  typical  of  turbulent  attached  flow  with  a  boundary  layer  height  of  5%  chord.  With  control  on,  the  existence 
of  the  wall-jet  yields  a  peak  velocity  exceeding  2Uoa  and  a  wall  velocity  gradient  higher  than  the  no-control  turbulent 
flow.  The  peak  is  slightly  larger  for  Re  =  45,  000  and  for  both  cases  u  relaxes  to  ux  =  L'^cost  1 5°)  in  the  far  field. 
Comparison  of  surface  Cp  at  Re  =  90, 000  with  and  without  control,  Fig.l4d,  indicates  that  the  effect  of  the  actuator 
is  relatively  small,  with  a  modest  reduction  of  pressure  over  much  of  the  upper  surface  and  a  small  spike  near  the  point 
of  application  of  the  actuator.  For  this  actuator  location,  the  main  impact  of  control  is  to  replace  a  turbulent  attached 
boundary  layer  with  a  laminar  attached  wall-jet. 


The  two  prominent  fluid  dynamic  effects,  turbulence  and  separation,  can  be  coupled  to  actuator  positioning  and 
orientation  to  yield  several  different  phenomena.  In  Ref.  32,  a  “backward”-oriented  actuator,  i.e.,  where  the  actuator  is 
placed  so  that  the  embedded  electrode  is  upstream  of  the  exposed  electrode  and  the  body  force  points  predominantly 
against  the  flow,  is  employed  to  induce  attachment  by  triggering  turbulence  in  a  separated  flow.  In  an  analogous 
fashion,  the  “forward”  facing  actuator  employed  above  may  be  deployed  to  trigger  separation  by  suitable  placement 
in  a  turbulent  flow  by  laminarizing  parts  of  the  flow  field.  To  demonstrate  this,  the  actuator  in  the  previous  case  was 
moved  from  near  the  leading  edge  to  the  point  of  maximum  thickness  in  the  turbulent  attached  flow,  as  depicted  in 
Fig.  15a.  Instantaneous  velocity  and  vorticity  magnitude  contours  are  plotted  in  Fig.  15b  and  c  respectively  after  an 
asymptotic  state  is  reached.  Comparison  with  the  baseline  flow  presented  earlier  in  Fig.  12  indicates  that  the  actuator 
has  triggered  separation  near  the  leading  edge:  in  fact  gross  manifestations  of  separation  appear  shortly  after  the 
actuator  is  turned  on.  The  flow  field  takes  on  the  character  of  the  no-control  lower  Reynolds  case  in  other  aspects 
as  well:  the  instantaneous  vorticity  contours  of  Fig.  15c  display  a  trailing  edge  vortical  structure  that  resembles  that 
observed  at  Re  =  45,  000  (see  Ref.  23).  These  results  are  significant  in  several  respects.  They  demonstrate  the 
potential  for  global  actuator  impact  and  highlight  the  need  to  consider  off-design  performance.  The  dynamics  indicate 
further  that  inviscid  or  2-D  simulation  procedures  are  unsuitable  for  such  studies. 

4.  Control  of  separated  flow  at  high  Reynolds  number,  Re  =  90, 000,  aa  =  18° 

As  noted  earlier,  increasing  the  Reynolds  number  from  45,  000  to  90, 000  at  an  angle  of  attack  of  15°  results  in  a 
turbulent  attached  flow.  In  order  to  examine  the  effect  of  actuation  in  a  separated  higher  Reynolds  number  flow,  the 
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a)  Force  location 


b)  Instantaneous  u 


c)  Instantaneous  |cj| 


d)  Mean  u 


Figure  15.  Effect  of  actuator  located  at  maximum  thickness  location  at  Re  =  90, 000,  aa  =  15° 


angle  of  attack  was  increased  to  18°.  The  main  features  of  the  flow  without  control  are  shown  in  Fig.  16.  The  effect  of 
angle-of-attack  may  be  discerned  by  comparison  with  Fig.  12,  which  shows  results  at  the  same  Reynolds  number  and 
aa  =  15°.  The  instantaneous  zt-velocity  and  vorticity  indicate  that  the  flow  is  turbulent  over  the  upper  surface  of  the 
wing  section,  while  the  time-averaged  velocity.  Fig.  16c  shows  that  separation  occurs  at  about  x/c  ~  0.5. 


a)  Instantaneous  u 


b)  Instantaneous  |ca| 

Figure  16.  Flow  features  at  Re  =  90,  000,  aa  =  18°  (no  control) 


c)  Mean  u 


The  effect  of  control  on  this  separated  flow,  with  the  actuator  placed  near  the  leading  edge,  is  shown  in  Fig.  17. 
In  the  initial  transients  (not  shown),  the  wall  jet  is  established  over  the  entire  wing  surface.  This  “impulse”  effect, 
discussed  in  Ref.  23,  is  short  lived  however,  as  instabilities  arise  in  the  shear  layer  downstream  of  the  trailing  edge. 
These  disturbances  subsequently  propagate  upstream  and  yield  a  turbulent  boundary  layer  over  the  wing  as  shown 
with  instantaneous  zt-velocity  contours  in  Fig.  17a  .  The  mean  zz-velocity  contours.  Fig.  17b  show  the  wall  jet  as  a 
clear  feature  until  X  ~  0.8,  beyond  which  turbulent  diffusion  yields  a  thicker  disturbed  region  but  with  lower  peak 
velocities.  Figure  17c  exhibits  the  mean  velocity  profiles  at  the  mid-chord  for  three  cases  in  which  the  Reynolds 
number  and  angle-of-attack  are  varied  independently  but  with  the  same  coefficient  Dc  =  2400.  All  cases  exhibit 
the  wall-jet  structure,  and,  since  Dc  scales  the  electric  with  the  inertial  force,  peak  velocities  in  the  jet  are  similar, 
reaching  about  2.5£/oo.  The  principal  differences  are  observed  in  the  near-wall  profile  fullness  which  is  highest  in  the 
high-Reynolds  number  and  angle-of-attack  case  and  the  wall-jet  remains  turbulent. 


c)  //.-pro (lies  at  mid-chord 

Figure  17.  Effect  of  control  with  actuator  located  near  leading  edge  at  aa  =  18°  and  Re  =  90,  000 
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5.  Effect  of  actuator  strength 

The  effect  of  actuator  strength  on  the  flow  field  is  examined  for  Re  =  90, 000  and  angle-of-attack  of  18°  by  reducing 
Dc  from  2400  for  the  previous  cases  to  half  that.  Results  shown  in  Fig.  1 8  indicate  that  there  is  no  significant  qualitative 
impact  associated  with  the  reduction.  The  instantaneous  //-velocity  shows  that  the  boundary  layer  remains  turbulent, 
similar  to  the  situation  depicted  in  Fig.  17a.  Further,  the  mean  velocity.  Fig.  18b  indicates  that  the  flow  remains  attached 
though  the  prominence  of  the  wall-jet  is  clearly  diminished  when  Dc  is  lowered.  Although  not  shown,  each  turbulent 
attached  flow  exhibits  similar  surface  Cp  profiles.  The  quantitative  differences  are  displayed  in  the  velocity  profiles 
at  mid-chord,  shown  in  Fig.  18c.  The  lower  forcing  yields  a  wall-jet  with  a  significantly  smaller  peak  velocity  l.bUoo 
versus  2.5£/oc.  These  results  suggest  that  the  effect  of  turbulence  in  establishing  attachment  augments  that  of  the 
development  of  a  wall-jet. 


a)  Instantaneous  //-velocity  b)  Mean  //.-velocity  c)  Comparison  of  //-profiles 

Figure  18.  Effect  of  actuator  strength  at  angle-of-attack  18°  and  Re  =  90,  000 


C.  Effect  of  Duty  Cycle 

Duty  cycles,  in  which  the  actuator  is  switched  on  intermittently  in  a  usually  regular  pattern,  have  been  explored  to 
enhance  effectiveness  in  the  context  of  dielectric  barrier  discharges  either  to  reduce  power  requirement  or  to  influence 
fluid  dynamic  instabilities  to  accomplish  specific  goals.10  For  the  wing  section,  both  objectives  can  be  achieved 
through  the  combined  effects  of  optimizing  near  wall  momentum  injection  and  triggering  of  transition  to  turbulence. 
To  examine  these  phenomena,  the  phenomenological  force  model  was  employed  for  the  stalled  Re  =  45, 000,  0  =  15° 
condition  with  two  different  inter-pulse  periods,  Tp  =  0.7 Tc  and  Tp  =  0.25TC,  of  same  duty  cycle  (20%)  and  power 
requirement.  The  first  value  0.7TC  represents  the  time  required  for  perceptible  manifestation  of  separation  to  appear 


Figure  19.  Amplitude  variation  in  duty  cycle,  Tp  =  0.7 Tc 


after  cessation  of  the  body  force  in  a  fully  attached  situation  (see  Ref.  23  for  transient  effects  after  the  actuator  is 
switched  off).  Results  obtained  with  this  forcing  suggest  exploration  with  the  lower  value.  The  amplitude  of  the 
imposed  force  for  the  Tp  =  0.7 Tc  case  is  shown  in  Fig.  19  as  a  function  of  time.  With  duty  cycle  of  20%,  the  actuator 
is  on  for  0.14TC  in  each  cycle.  The  start  and  end  of  the  on-period  is  smoothed  with  a  fifth  order  polynomial: 

A(t)  =  at5  +  btA  +  ct3  +  dt2  +  et  +  f  (12) 
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For  the  ramp-up  segment,  say  between  t\  and  i2,  the  coefficients  are  obtained  by  enforcing  A  ( 1 1 )  =  0,  A(t2)  =  1  and 
dA/dt  =  0  and  d2A/dt2  =  0  at  t±  and  i2: 


a  =  — 


D 


d  = 


_  30tit2(ti+t2) 

D 


b  = 
e  = 


15(ti+t2) 

D 

30 t*t% 
D 


-10(t?+4tit2+t|) 

D 

D 


D  —  (ti  —  f2)5 


(13) 


10%  of  the  duty  cycle  is  employed  for  ramp-up  and  ramp-down  (i.e.,  0.014 Tc  for  Tp  =  0.77)  and  0.005  for  Tp  = 
0.25  Tc). 

The  response  may  be  naturally  divided  into  two  time  ranges:  the  initial  transients  and  the  asymptotic  state.  The 
former  are  shown  for  Tp  =  0.7T)  in  Fig.  20  which  depicts  instantaneous  u  contours  at  selected  points  in  the  first  cycle, 
with  time  being  measured  from  the  onset  of  control.  The  initial  condition,  Fig.  20a,  corresponds  to  an  instantaneous 
pattern  of  the  separated  flow  with  a  shear  layer  that  rapidly  becomes  turbulent.  The  onset  of  the  body  force  interrupts 
the  formation  of  the  shear  layer,  together  with  the  eruption  of  a  sequence  of  vortical  structures,  similar  to  those 
described  in  Ref.  32.  By  the  time  instant  of  Fig.  20c,  the  force  has  been  switched  off.  Subsequently,  the  initial 
structures  in  the  disturbed  region  are  essentially  convected  downstream  without  the  sustaining  effect  of  the  vortex 
sheet.  At  the  end  of  the  interpulse  period,  incipient  separation  is  evident  near  the  leading  edge,  but  its  development  is 
inhibited  by  the  onset  of  the  next  actuation  cycle. 


Figure  20.  Initial  transients  after  actuation  with  20%  duty  cycle  and  interpulse  period  0.7 Tc.  Time  measured  after  onset  of  actuation. 


The  sequence  of  events  when  the  asymptotic  state  is  reached  is  shown  in  Fig.21.  Figure  21a  depicts  the  situation 
at  roughly  the  mid-point  of  the  off  segment  of  the  cycle.  A  major  feature  in  each  figure  is  the  existence  of  relatively 
large  coherent  structures  outside  of  the  near  wall  jet  region  which  are  convected.  This  is  observed  to  be  the  direct 
result  of  the  impulsive  nature  of  the  forcing.  Upon  switching  on  the  actuator  a  set  of  vortices  is  initiated  which  lift 
from  the  surface  of  the  wing.  In  the  figures,  two  such  structures  arising  during  successive  on-phases  of  the  actuator 
have  been  identified  as  A  and  B  respectively.  These  structures  diffuse  as  they  are  convected  downstream,  partly  as  a 
consequence  of  the  coarser  mesh,  but  maintain  their  overall  form  and  can  be  identified  even  near  the  trailing  edge.  By 
frame  (e),  separation  inception  is  apparent  near  the  leading  edge,  becoming  prominent  by  frame  (f).  In  frames  (g)  and 
(h),  the  actuator  is  on,  injecting  momentum  near  the  surface  and  initiating  the  formation  of  a  new  local  jet  preceeded 
by  the  vortex  pair  B,  which  follows  the  same  development  as  A,  thus  repeating  the  sequence. 

The  development  of  incipient  separation  in  the  previous  simulations  suggests  that  the  interpulse  period  of  0.7TC  be 
reduced.  To  explore  the  effect  of  shorter  bursts  of  more  frequent  actuation,  a  simulation  was  conducted  at  Tp  =  0.25 Tc. 
With  the  same  duty  cycle  and  ramp-up  and  ramp-down  times  (as  a  ratio  of  the  duty  cycle)  the  power  requirements 
remain  the  same.  Results  are  shown  in  Fig. 22  after  a  periodic  asymptote  is  reached.  Note  that  the  time  elapsed 
between  frames  is  now  considerably  shorter  compared  to  that  in  Fig.  21.  The  flow  downstream  and  near  the  surface 
is  more  disturbed  than  in  the  longer  cycle.  The  results  show  that  the  eruption  of  vortical  structures  persists  in  this 
case.  Close  scrutiny  indicates  that  up  to  three  distinct  vortex  pairs  of  varying  strength  -  the  leading  pair  being  the 
most  distinct  -  are  generated  with  each  actuator  onset  event.  Overall,  the  size  of  these  structures  is  smaller  than  in  the 
previous  case  and  their  trajectory  remains  relatively  closer  to  the  surface. 
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e)  T  =  10.475TC 


f)  T  =  10.5625TC  g)  T  =  10.65TC 

Figure  21.  Asymptotic  response  with  20%  duty  cycle  of  interpulse  period  of  0.7 Tc 


h)T  =  10.7375TC 


e)T=  10.15TC 


f )T  =  10.187571 


g)  T  =  10.225 Tc 


Figure  22.  Asymptotic  response  with  20%  duty  cycle  of  interpulse  period  of  0.2571 


h)  T  =  10.2625TC 

:  u-velocity 


The  time-averaged  response  to  actuation  with  the  two  interpulse  periods  is  shown  in  Fig.  23.  Figures  23a  and 
b,  depict  the  u-velocity  after  an  asymptotic  state  is  reached.  The  development  of  the  initial  stages  of  separation  for 
the  longer  interpulse  period  case  is  reflected  in  a  bulging  of  the  boundary  layer  near  the  leading  edge  (see  inset)  and 
the  boundary  layer  displays  rapid  growth  aft  of  the  mid-chord  region  because  of  turbulence.  The  mean  flow  remains 
attached  over  the  upper  surface.  In  contrast,  the  smaller  interpulse  period  case  exhibits  separation  in  the  mid-chord 
region  but  the  smaller  boundary  layer  height  near  the  trailing  edge  reflects  the  closer  proximity  of  the  vortical  structures 
to  the  wing  surface.  Mean  velocity  profiles  are  compared  at  mid-chord  and  one  chord-length  downstream  of  the  trailing 
edge  in  Fig.  23c.  At  mid-chord,  the  shorter  interpulse  period  yields  a  shallower  velocity  profile  but  a  sharper  and  higher 
peak  ( u/Uoo  ~  1.4),  while  the  larger  period  case  exhibits  a  more  diffuse  distribution.  However,  these  differences  are 
mitigated  downstream,  where  similar  velocity  defects  are  observed  in  the  wake.  Mean  Cp  profiles  plotted  in  Fig.  23d 
indicate  that  with  the  exception  of  the  region  in  the  immediate  vicinity  of  the  actuator,  the  mean  pressure  fields  are 
similar  with  the  two  different  interpulse  periods. 
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a)  Mean  u-velocity  for  Tp  = 
0.7  Tc 


v  m 


b)  Mean  v/.- velocity  for  Tp  =  c)  Mean  u  profiles 

0.25TC 

Figure  23.  Effect  of  interpulse  period  variation 

V.  Conclusions 


Full  Navier-Stokes  simulations  are  reported  on  plasma-based  flow  control  of  a  stalled  NACA  0015  wing  section 
at  Reynolds  numbers  of  45, 000  and  90,  000  and  angles  of  attack  15°  and  18°.  The  time-accurate  force  field  obtained 
from  a  self-consistent  simulation  at  5 kHz  and  2 kV  rms  is  suitably  scaled  and  employed  to  determine  the  response 
of  the  flow.  The  unsteady  asymptotic  pattern  as  well  as  the  mean  flow  are  described  in  terms  of  the  small-scale 
structures,  acoustic  signals  and  the  development  of  a  wall-jet.  Separation  location  moves  downstream  with  perceptible 
detachment  occurring  at  80%  chord.  In  many  key  respects,  the  time-averaged  flowfield  is  similar  to  that  obtained  from 
a  steady  time-averaged  force  and  that  reported  with  a  phenomenological  model  in  an  earlier  publication,  though  more 
investigation  is  warranted. 

This  simpler  averaged  model  is  employed  with  3-D  direct  numerical  simulations  to  deduce  Reynolds  number,  force 
scaling  and  duty  cycle  effects.  At  higher  Reynolds  number,  90, 000,  the  no-control  flow  at  15°  angle-of-attack  becomes 
turbulent  immediately  downstream  of  the  leading  edge  and  remains  attached  in  the  mean.  The  actuator  laminarizes  the 
flow  with  the  formation  of  a  wall-jet.  When  the  actuator  location  is  moved  to  the  point  of  maximum  airfoil  thickness, 
the  flow  near  the  leading  edge  becomes  laminar  and  separates.  To  explore  stall  control  at  Re  =  90,  000,  the  angle  of 
attack  was  raised  to  18°.  Actuation  at  the  leading  edge  recovers  the  attached  laminar  wall  jet  on  the  upper  surface. 
The  effect  of  actuator  force  strength  on  the  response  is  examined  and  the  combined  effects  of  the  laminar  wall  jet 
and  subsequent  turbulence  breakdown  are  described.  Finally,  the  effect  of  pulsing  is  documented  at  Reynolds  number 
45,  000  with  20%  duty  cycle  and  interpulse  periods  of  0.7  and  0.25  characteristic  times  respectively. 
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